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Abstract 

The influence of cold and heated secondary 
flow on the instability of a two-stream, coplanar 
jet having a 0.7 Mach number heated primary jet 
for a nominal fan to primary velocity ratio of 
0.68 was investigated by means of inviscid lin- 
earized stability theory. The instability prop- 
erties of spatially growing axisymmetric and first 
order azimuthal disturbances were studied. It was 
found that the instability characteristics of the 
two-stream jet with a velocity ratio of 0.68 are 
very different from those of a single stream jet, 
and a two-stream, coplanar jet having a 0.9 Mach 
number heated primary jet and a cold secondary jet 
for a fan to primary velocity ratio of 0.30. For 
X/D = 1 and in comparison to the case where the 
velocity ratio was 0.3, the presence of the fan 
stream with a velocity ratio of 0.68 enhanced the 
instability of the jet and increased the unstable 
frequency range. However, the axisymmetric mode 
(m = 0) and the first order azimuthal mode (m = 1) 
have similar spatial growth rates where the veloc- 
ity ratio is 0.68 while for a velocity ratio of 
0.3 the growth rate of the first order azimuthal 
mode (m = 1) is greater. Comparing the cold and 
hot secondary flow results showed that for a 
velocity ratio of 0.68 the growth rate is greater 
for the case with cold secondary flow. 

Nomenclature 


a speed of sound, m/s 

B see Eq. (1) 

Cph phase speed, m/s 

F(r,B,s) radial profile term 
f eigenfunction 

G(r) velocity profile shape 

H(a ) eigenvalue equation 

i V* 

m azimuthal wavenumber 

p(r) pressure 

R radius 

r radial position 

T(r) static temperature 

t time, sec 

U c ix mean velocity on centerline at a given 

x/ Dp 

u,v,w velocity components 


x axial position 

a complex eigenvalue wavenumber 

3 radian frequency 

au velocity difference, U c ] x - !)«> 

5 see Eq. (1) 

e momentum boundary layer thickness 

x see Eq. (8) 

p density 

cylindrical angular coordinate 
Subscripts : 

c compressible 

clx center line at a given x/D p position 

e exterior solution 

h meeting point 

i imaginary 

p primary jet 

r real 

s secondary or fan jet 

t temperature related 

0 interior solution 

® ambient 

Superscript : 

( )' small disturbance 

(-) mean value 

Introduction 

For many future aircraft, improvements in 
operation will be in part due to control and modi- 
fication of jet exhaust velocity and/or temperature 
characteristics. For example, the jet velocity and 
temperature impinging in the flap surfaces of a 
short takeoff and landing aircraft employing an 
under-the-wing blown flap for lift augmentation dur- 
ing takeoff and landing must be controlled in order 
to maintain flap loads and surface temperatures 
within acceptable limits for reasons of structural 
integrity and high-lift performance. Commercial 
supersonic transports and transatmospheric vehicles 
are other examples of future aircraft for which 
changes in the jet exhaust velocity and temperature 
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may be important since the jet exhaust noise must 
be reduced to achieve acceptable community noise 
levels, Coannular nozzles, mechanical suppressors, 
and thermal acoustic shields have successfully 
reduced jet noise. 


only the most rapidly growing spatial modes, 
information will be presented only on this type 
of disturbance. 

Theory 


Temperature and velocity control is usually 
achieved through jet mixing with the surrounding 
medium. Normally, mixing is achieved using complex 
nozzle or shroud/ejector geometries with consequent 
weight penalties. Enhanced mixing using acoustic 
or aero/mechanical perturbing devices may reduce 
this weight penalty. Other situations where acous- 
tic and aero/mechanical excitation may be benefi- 
cial are discussed in Ref. 1. Due to these 
possible benefits the NASA Lewis Research Center 
has been supporting experimental, numerical, and 
theoretical research programs in this area.^“° 

The effect of excitation on centerline veloc- 
ity and static temperature decay for a single plume 
jet was estimated using correlation methods in 
Ref. 9. Some of the data used in these correla- 
tions is from the results of the comprehensive 
experimental and theoretical research program 
presented in Refs. 2 to 6 which included both hot 
and cold flows. Additional cold flow results are 
reported in Ref. 8. As part of a program to obtain 
similar information for mixing jets, the unexcited 
velocity and temperature characteristics of two- 
stream coplanar jet exhaust plumes was investigated 
in Ref. 10. In addition to information on axial 
velocity and temperature decay. Ref. 10 also pre- 
sents typical radial velocity and temperature pro- 
files. In an earlier paper, the influence of the 
second stream on the instability of a two-stream, 
coplanar jet with only the primary stream heated 
for a nominal secondary to primary velocity ratio 
of 0.3 was investigated.^ In the present paper, 
the radial velocity and temperature profiles for a 
secondary to primary velocity ratio of 0.68 for hot 
and cold secondary flow are used to determine the 
phase speed and spatial growth rate of pressure 
instability waves. 

An excellent review of recent research on the 
instability of spatially growing disturbances in 
the free boundary layer of a jet having a single 
inflection point was presented by Michalke.* 2 
The two-stream coplanar jet studied herein has two 
inflection points. As far as the author can deter- 
mine, little work has been done to study free shear 
layers with two inflection points. ^ However, the 
instability of a circular incompressible jet with 
external flow was analyzed by means of linearized 
stability theory by Michalke and Herman.^ For 
the Mach number range of the present study (around 
M = 0.7) the flow is compressible. The instability 
of spatially growing disturbances propagating isen- 
tropically in a circular jet for flow in this Mach 
number range was studied by Michalke^ and 

Morris. 16 y 

In this study viscous effects and effects due 
to the slowly diverging jet flow!' - * 9 are 
neglected. Furthermore, the analysis follows that 
of Michalkel 5 where it was assumed that disturb- 
ances propagate isentropically. Because the large 
scale structure of turbulence in a circular jet is 
dominated by the axisymmetric and first order azi- 
muthal components, ^6 the stability calculations 
are restricted to these two components. Since the 
long term behavior of the disturbances depends on 


Jet Velocity and Temperature Characteristics 


In order to determine the stability of the 
disturbances, the radial distribution of mean 
velocity and temperature in the two-stream axisym- 
metric jet must be specified. The conditions 
selected for study are among the cases presented 
in Ref. 10 as being typical. For these cases the 
nominal primary jet Mach number is 0.7, the nomi- 
nal fan-to-primary velocity ratio is 0.68 and the 
primary stream is heated. The nozzle used for 
this case had a secondary to primary flow area 
ratio of 1.9 since the primary nozzle diameter is 
10 cm, the secondary nozzle diameter is 17.6 cm, 
and the primary nozzle thickness is 0.46 cm. The 
actual velocities and temperatures are shown in 
Table I. In Table 1(a) the case with the cold 
secondary is shown while in Table 1(b) the case 
with the hot secondary is presented. 


The velocity and temperature profiles were 
constructed using a basic profile shape. The 
same method was used in Ref. 11. The basic pro- 
file used is 


F(r,B,6 ) 


1 

2 


1 + tanh 



( 1 ) 


where the constants B and 6 are determined 
from experimental data. The constants used were 
determined by curve fitting data given in Ref. 10 
and they are shown in Table I. Measured and cal- 
culated velocity and temperature profiles are 
compared in Figs. 1 and 2. 


It was found that a reasonable curve fit to 
the experimental velocity profile could be 
achieved by adding together two separate profiles 
specified by Eq. (1) with appropriate constants 
and appropriate weighting coefficients. The 
resulting equation for the curve fit to the veloc- 
ity profile is 

G u (r) = “ uj) F ( r>B u p ’ s u p ) + uj ( F r ’Vu s ) (2) 

Near the nozzle where the velocity profile for the 
primary and secondary jets have a top hat profile 
the values of B are large (B.. = 10 and 

P 

By = 20) and Gy(0) is unity. However, at large 

x/D the profiles are broad and the values of B 
are nearer unity. In this region at r = 0 the 
value of G ( 0) is not unity. Since the actual 
velocity profile must at r = 0 have a value of 
unity, the basic shape is normalized by dividing 
by Gy ( 0) . Thus the actual velocity profile 
used is 


u (r) Gj(r) 



where U c i x is the mean velocity on the center- The disturbance is assumed to be a spatially 

line at a given x/D p position. growing wave of the type 


The equation used to curve fit the experi- 
mental temperature profile for the case with a 
cold secondary is 

Mr) = p (r,B ,« T ) (4) 

1 1 P 1 P 


The constants used were determined by curve fit- 
ting data given in Ref. 10 and they are shown in 
Table I. The temperature profile used for the 
case with a hot secondary is similar to that used 
for the velocity profile 



(5) 


Again, just as for the velocity profiles the 
actual velocity profile used for this case is 
normalized by Gj(0) . 


T(r) - T^ Mr) 


( 6 ) 


Again the constants used were determined by curve 
fitting data given in Ref. 10 and they are shown 
i n Table I . 


[u\v , ,w\p , ,p'] = [u(r),v(r),w(r),p(r),p(r)] 

exp[i(ax - m<p - et)] (9) 

where the radian frequency, a, and the integer 
azimuthal wave number, m, are real, while the 
eigenvalue wavenumber, o, is complex. The real 
value of a, a r is the axial wave number. The 
imaginary value of a determines the spatial 
growth rate, -a-j. If -a-; is greater than zero 
the wave defined by Eq. (9) grows as it propagates. 
The propagation phase velocity c p h is given by 

8 / op * 


Using Eq. (9), the disturbance equations can 
be reduced to a single equation for the pressure 
perturbation p(r) : 


L(p) 



where 


2 ) dU(r) dp 

[U(D - B/a] ) Ar Ar 
[a 2 x 2 (r) + m 2 /r 2 ] p = 0 


( 10 ) 


x 2 (r) = 1 - 



2 


The momentum loss boundary layer thickness, 
e c , is used as the characteristic of the free jet 
shear layer. For compressible flow, it is defined 
as follows: 



Note, however, that the incompressible momentum 
boundary layer thickness defined as 


and a is the acoustic speed of sound. Equations 
of this form for the fluctuating pressure are well 
known and given in Refs. 15 and 16. 

Near the jet axis (r -► 0) and in the ambient 
fluid far from the jet (r -► ®) the velocity profile, 
U(r), is constant and dU(r)/dr is zero. In 
these regions Eq. (10) reduces to a form of the 
Bessel equation. Consequently, the asymptotic 
solutions to Eq. 10 are given by the modified 
Bessel functions I m and K m of order m. The 
boundary conditions for the pressure require that 
p(0) has to be bounded and p(«>) has to be zero. 
Hence, near the jet axis (r + 0) 



is more frequently used to characterize the thick- 
ness of the velocity shear layer for jet shear 
layer instability studies even for compressible 
flow. 12 Consequently, it is used again in this 
paper. 

Linearized disturbance equations 

The disturbance equations are derived by 
introducing small disturbances u ' , v‘ ,w' ,p' , p ' 
into the equations of momentum and mass conserva- 
tion and keeping only the linear terms. The 
equations are linearized about the basic jet flow 
which is assumed to be locally parallel with U(r) 
being the axial velocity component of the undis- 
turbed flow. The Reynolds number of theflow is 
assumed to be large so that the flow is inviscid. 
The flow is also assumed to be isentropic. 


p(r) s I m [ax(r)r] (11) 

and in the ambient fluid (r ♦ «>) 

p(r ) * K m (ax(r)r) (12) 

The numerical procedure used to solve the 
eigenvalue problem is a shooting method similar to 
that described in Ref. 15. Another method based 
on a finite difference scheme has been recently 
developed. Both methods give similar results 
when applied to the same test problems. However, 
near the neutral stability point, convergence is a 
problem with the shooting method. The step size 
must be reduced until one is sure that the result 
is independent of this parameter in order to be 
sure the shooting method has produced a valid 
result. 

The differential Eq. (10) is scaled using the 
incompressible momentum boundary shear layer thick- 
ness, e, and the velocity difference aU, between 
the jet centerline velocity at a particular x/D p , 
U c ] x , and the ambient jet velocity, U^. For the 
case studied, aU is the jet centerline velocity. 
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U c ix» since the ambient velocity is zero. The 
resulting differential equation then depends on two 
normalized parameters, the dimensionless spatial 
growth rate -a^e and the dimensionless frequency, 
se/Al). 

The complex differential equation written as a 
first order vector system is integrated numerically 
by a Runge-Kutta procedure. The infinite integra- 
tion region is divided into two finite regions, an 
inner region r 0 < r < and an outer region 
rh < r < r ro . The values of r 0 and r w are 
chosen such that dU/dr is small and the integra- 
tion of Eq. (10) can be started at these points. 

In the inner region 

p 0 ( p ) = f 0 ( r ) ( i3 ) 


where fo(r) represents the numerical Runge-Kutta 
solution in the inner region. In the inner region 
the numerical integration goes from r^ toward 
r^. Near the origin the proper initial condition 
for the Runge-Kutta integration from r 0 to r^ 
is from Eq. (11) 

Pq^O 5 = I m CaX(r 0 ,r 0 ] (14) 


In the outer region the numerical Runge-Kutta 
solution is given by 


p o (r o 


- C eV r > 


(15) 


In contrast to the inner region, this Runge-Kutta 
solution is obtained by integrating from r ro 
toward the origin. Far from the jet from Eq. (12) 
the proper initial condition for the Runge-Kutta 
integration from r„ to r^ is 

•VO = W“ x(r >J os) 


Results 


Calculated axisymmetric (m = 0) and azimutha' 
(m = 1) spatial growth rates and phase velocities 
at x/Dp = 1 for a two-stream coplanar jet with 
a secondary to primary velocity ratio of 0.68 and 
for both cold and hot secondary flows are shown in 
Figs. 3 and 4. Also shown are results presented 
in Ref. 11 obtained at x/D = 1 for m = 0 and 
m = 1 for a two-stream coplanar jet with a sec- 
ondary to primary velocity ratio 0.3 with cold 
secondary flow. At x/ Dp = 1 the boundary layer 
thickness for the case where the velocity ratio is 
0.3 is 0.899 cm. The boundary layer thickness for 
the cases where the velocity ratio is 0.68 is 
5.11 cm when the secondary flow is cold and 3.08 cm 
when it is hot. The nozzle geometry for all cases 
is identical. Figure 3 shows the dimensionless 
spatial growth rate, ~a-j,e as a function of the 
dimensionless frequency, bo/au, for the axisym- 
metric (m = 0) disturbance and for the first order 
azimuthal (m = 1) disturbance. At all frequencies, 
for a velocity ratio of 0.68, the axisymmetric and 
the first order azimuthal disturbances have nearly 
identical growth rates. Furthermore, the magnitude 
of the growth rate for the case of cold secondary 
flow is much greater than that of the case with 
hot secondary flow. Note that the spatial growth 
rate with a velocity ratio of 0.68 is greater than 
that for a velocity ratio of 0.3. 

For the cold secondary rase with a jet veloc- 
ity ratio of 0.3, the preferred excitation Strouhal 
number is 0.05. However, the growth rate is very 
low. For this case it is better to excite the 
first order azimuthal mode at a Strouhal number of 
0.45. For the case with a velocity ratio of 0.68 
and for both modes, the case with a heated second- 
ary flow has its largest growth rate at a Strouhal 
number of 0.2. However, the case with the cold 
secondary flow has its largest growth rate at a 
Strouhal number of 0.4. 


The eigenvalue is determined by an iterative 
procedure. The eigenvalue is the value of a 
which makes p and dp/dr continuous at r = r^. 
Hence, the eigenvalue equation for a is 


H(a) f 0 ( r h i 


df e (r b ) f e (r h )df 0 (r h) 


Kr r 


“3 F 


= 0 


(17) 


where fo and C e f e are the solutions in the 
inner and outer regions. The constant C e is 
determined by the condition that at r^ the inner 
and outer pressure eigenfunction solutions have 
the same value. 


The eigenvalue is determined in a two step 
procedure. The first step uses a technique 
employed in Ref. 22. The technique is based on a 
theorem in complex variable theory which states 
that the number of zeros within a closed contour 
equals the net multiples of 2tt by which the 
phase angle of H(a) changes around the contour. 23 
This method is extremely powerful since by using 
it, one can be sure one has found the most unstable 
eigenvalue. The method also yields an accurate 
initial guess. The initial guess is then used by 
a least squares function minimization iterative 
procedure, which finds the eigenvalue. 


In Fig. 4 the axial phase velocity, 

Cph = 8/a r , for X/D =1 is plotted as a function 
of the frequency for m = 0 and m = 1. For a 
single jet the phase velocity of the axisymmetric 
disturbance (m = 0) always decreases with __ 
increasing frequency from the jet velocity U c ] x 
at zero frequency. Also, for a single jet the 
phase velocity of the first azimuthal disturbance 
is always less than the phase velocity of the 
axisymmetric one. 6 For the two-stream axisym- 
metric jet, the frequency dependence of the phase 
velocities of the axisymmetric and first azimuthal 
disturbances are similar to those of a single jet. 6 
For a velocity ratio of 0.68, the phase velocity 
is lower if the secondary is hot for the axisym- 
metric and first order azimuthal case. Note that 
the phase velocity for the case with a velocity 
ratio of 0.68 is lower than for the case where 
the velocity ratio is 0.3. 

Concluding Remarks 

It has been shown that the instability char- 
acteristics of the two-stream axisymmetric jet 
with a secondary to primary velocity ratio of 0.68 
are different from those obtained for a secondary 
to primary velocity ratio of 0.3. For example, 
for cold secondary flow with a jet velocity ratio 
of 0.3, the growth rate of the axisymmetric mode 
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(m = 0) is much less than the growth rate of the 
first order azimuthal mode (m = 1). While for 
cold or hot secondary flow with a jet velocity 
ratio of 0.68. The growth rates of the axisym- 
metric mode (m = 0) and first order azimuthal 
mode (m = 1) are equal . 

The presence of the heated secondary stream 
reduced the instability of both modes. For a two- 
stream, coplanar jet operating at a secondary to 
primary velocity ratio of 0.68 aero/mechanical 
excitation can be designed to excite either mode. 
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TABLE I. - JET TEST CONDITIONS AND PROFILE PARAMETERS 
(a) Cold Secondary Flow 


Exit Mach number 0.69 

Primary jet nozzle diameter, cm 10 

Secondary or fan jet nozzle diameter, cm . . . 17.6 

Primary jet nozzle wall thickness, cm 0.46 

Secondary to primary nozzle area ratio .... 1.9 

Primary jet exit velocity, m/s 459 

Secondary jet exit velocity, m/s 314 

Primary jet exit temperature, K 1092 

Secondary jet exit temperature, K 248 


X/Dp 

Rp/e 

e 

cm. 

e c, 

cm 

D clx, 

m/s 

T clx, 

K 

B U 

Up 

«n 

Up 

\ 

\ 

b t 

'p 

Up 

1 

1 

5.11 

1.32 

459 

1092 

10 

0.8 

20 

1.9 

6 

0.7 

2 




459 

1081 

5 

.6 

8 

1.8 

6 

.8 

4 




433 

1062 

5 

.6 

5 

2 

3 

.5 

8 




369 

693 

2 

.6 

3 

2.2 

.8 

.3 

12.75 




290 

403 

1 

.7 

4 

.3 

.8 

.3 


(b) Hot Secondary Flow 


Exit Mach number 0.71 

Primary jet nozzle diameter, cm 10 

Secondary or fan jet nozzle diameter, cm . . . 17.6 

Primary jet nozzle wall thickness, cm 0.46 

Secondary to primary nozzle area ratio .... 1.9 

Primary jet exit velocity, m/s 459 

Secondary jet exit velocity, m/s ...... . 314 

Primary jet exit temperature, K 1050 

Secondary jet exit temperature, K 545 


X/Dp 

Rp/e 

6 

cm. 

6 C, 

cm 

D clx, 

m/s 

T clx, 

K 

\ 

\ 

\ 


b t 

'P 

«t 

'P 

% 

\ 

1 

1 

3.08 

1.19 

459 

1050 

20 

0.9 

10 

1.8 

10 

0.9 

10 

1.7 

2 




459 

1050 

5 

1.1 

8 

1.9 

5 

.9 

5 

1.7 

4 




459 

1050 

5 

.9 

2 

1.8 

2 

.7 

4 

1.7 

8 




415 

883 

2 

.6 

2 

1.6 

2 

.7 

2 

2 

12.75 




274 

633 

3 

1.3 

1 

2.3 

i 

.7 

1 

.7 

2 










pe/Au 

(B)m= i. 

FIGURE 3.- SPATIAL GROWTH RATE OF THE AX I SYMMETRIC (m = 0) 
AND FIRST AZIMUTHAL (m = 1) DISTURBANCE AT X/D = 1. 
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0.68 are very different from those of a single stream Jet, and a two-stream, 
coplanar Jet having a 0.9 Mach number heated primary Jet and a cold secondary Jet 
for a fan to primary velocity ratio of 0.30. For X/D = 1 and In comparison to 
the case where the velocity ratio was 0.3, the presence of the fan stream with a 
velocity ratio of 0.68 enhanced the Instability of the Jet and Increased the 
unstable frequency range. However, the axlsymmetrlc mode (m = 0) and the first 
order azimuthal mode (m = 1) have similar spatial growth rates where the velocity 
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order azimuthal mode (m = 1) Is greater. Comparing the cold and hot secondary 
flow results showed that for a velocity ratio of 0.68 the growth rate Is greater 
for the case with cold secondary flow. 
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